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Abstract. Violation of the axial symmetry of a magnetic 
field essentially modifies the physics of the plasma outflow in 
the magnetosphere of rotating objects. In comparison to the 
cixisymmetric outfiow, two new effects appear: more efficient 
magnetocentrifugal acceleration of the plasma along particu- 
lar field lines and generation of MHD waves. Here, we use an 
ideal MHD approximation to study the dynamics of a cold 
wind in the nonaxisymmetric magnetosphere. We obtain a self- 
consistent analytical solution of the problem of cold plasma 
outflow from a slowly rotating star with a slightly nonaxisym- 
metric magnetic field using perturbation theory. In the axisym- 
metric (monopole-like) magnetic field, the first term in the ex- 
pansion of the terminating energy of the plasma in powers of 
S7 is proportional to fi'', where fi is the angular velocity of the 
central source. Violation of the axial symmetry of the mag- 
netic field crucially changes this dependence. The first correc- 
tion to the energy of the plasma becomes proportional to Q. 
Efficient magnetocentrifugal acceleration occurs along the field 
lines curved initially in the direction of the rotation. I argue 
that all the necessary conditions for the efficient magnetocen- 
trifugal acceleration of the plasma exist in the radio pulsar 
magnetosphere. We calculated the first correction of the ro- 
tational losses due to the generation of the MHD waves and 
analysed the plasma acceleration by these waves. 
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1. Introduction 

The outflow of plasma from rotating magnetised objects is a 
widespread phenomenon in the Universe. It is observed in a 
wide range of astrophysical objects of different natures and 
scales. Some of these objects eject relativistic plasma. The 
Lorentz factor of plasma ejected by mi croqu asars and AGN 
is of the order 10 (Mirabel & Rodrigues 1996| , Begelman et al. 
1984, Pelletier et al. 1996), although there is evidence that 
in Blazars the plasma i s acc elerated to the Lorentz factor 
- 10" (Aharonian et al. |l999| ). Radio pulsars also accelerate 
the p lasma to the L orentz facto r ~ lO'' (Kennel & Coroniti 



1984, Coroniti 



1990 



Arons 



1996 ) . The problem of the acceler- 



ation of the plasma to such high energies remains one of the 
most important unsolved problems in modern relativistic as- 
trophysics. 



The acceleration of plasma can occur in principle due to 
different dissipative as well as nondissipative processes. Hy- 
pothesis have been proposed about the singular behaviour of 
the relativistic plasma outflow at the "l ight s urface" (1983) or 
at the light cylinder (Mestel & Shibata 1994) of radio pulsars. 
These dissipative structures could provide the plasma accelera- 
tion. However, recent direct numerical calculations of the rela- 



tivistic outflows in MHD approximation (B ogova lov 1997) and 
force- free approximation (Contpoulos et al. 1999) have demon- 



strated that the relativistic plasma flow is actually regular ev- 
erywhere. 

Another mechanism of acceleration due to dissipative pro- 
cesses was proposed by Michel (1982) and Coroniti (1990) in 
application to radio pulsars. Michel was the flrst who noticed 
that the number of charge carriers in the wind from radio pul- 
sars is not enough to sustain the necessary return current in 
the wind with a striped magnetic field formed at the outflow 
from the oblique rotator. Coroniti (1990) considered the re- 
connection process in this magnetic field and came to the con- 
clusion that it can provide the necessary acceleration of the 
winds from radio pulsars. Up to now this process was consid- 
ered as the only possible mechanism explaining pulsar wind 
acceleration. However, recently Luybarskii & Kirk (2000) have 
revisited the process of wind acceleration due to magnetic field 
reconnection. They concluded that well before the winds from 
pulsars such as the Crab pulsar reach high Lorentz factors they 
will be terminated by the interstellar medium. 

A lot of efforts has been spent to obtain the plasma accel- 
eration in frameworks of ideal MHD. There is a general opinion 
that in all astrophysical objects ejecting plasma, the magnetic 
field plays an essential role in the dynamics of the plasma. In 
particular, the magnetic field of a rotating object itself can pro- 
vide the acceleration of the plasma due to a m agnet ocentrifugal 
(below MC) mechanism (Blandford & Payne |l982| ). The ideal 
MHD outflows of plasma in different models and different ap- 
proximations were investiga ted in a series of w orks (Michel 
196E| , Sulk anen & Lovelace |l990| , Ferreir a |l997| , Vlahakis & 
Tsinganos 1998| , Bogovalov fc Ts inganos 199!: , Krasnopolsky 



et al. 1999, Ustyugova et al. 199t). It follows from these stud- 



ies that nonrelativistic plasma can be accelerated rather effi- 
ciently due to this mechanism. As far as relativistic plasma 
is concerned, the efficiency of the acceleration of this plasma 
appears insufficient for astrophysical applications. There were 
hopes that divergence of the magnetic fiux tubes somewhere 
beyond the light cylinder will result in more efficient accelera- 



2 



tion of the plasma due to the pre ssure gradient of the toroidal 
magn etic field (Begelman & Li 



1994, Takahashi & Shibata 



1995 ). The only known mechanism for this magnetic flux tube 



divergence is magnetic self-c oUimation of t he ma gnetised winds 
(Heyvaerts & Norman [l989|, Chiueh et al. |l99l|). This mecha- 



arman 1989, Chiueh et al. 
nism can, in principle, provide divergence of the field lines at 
the equatorial plane such that the poloidal field decreases with 
distance faster than r~^. This divergence can take place in the 
limiting range of distances from the source, since finally the 
wind expa nds ra dially near th e equ ator at r oo (Heyvaerts 
& Norman 



198£ , Chiueh et al. 



1991) 



Numerical self-consistent solution of the problem of steady- 
state plasma outflow in a model with the an initially monopole- 
like magnetic field shows that coUimation really increases the 
efficiency of nonrelativisti c plas ma acceleration at the equator 
(Bogovalov & Tsinganos 1999). However, coUimation of the 
relativistic plasma outfiow appe ars negligib lv srn all at large 
Lorentz factors (Bogovalov 1997 , Bogovalov 2001 ). Therefore, 
the relativistic plasma is practically not accelerated by the 
pressure gradient of the toroidal magnetic field. This conclu- 
sion is valid not only for the monopole-like model. 

The acceleration of the wind due to the pressure gradient 
of the toroidal magnetic field occurs at the large distances com- 
pared to the dimension of t he cen tral source (Begelman & Li 
1994 , Takahashi & Shibata 1998| ) . The stationary wind from 
any source is monopole-like at these distances. Therefore, the 
poloidal magnetic field also becomes monopole-like, since the 
magnetic field is frozen into the plasma. Thus, the solutions 
obtained in the monopole-like model actually describe general 
properties of any axisymmetric wind at large distances from 
the source. 

The application of the monopole-like model is not limited 
by large distances. This model also gives upper limits on the 
efficiency of the MC acceleration of the plasma in any other 
axisymmetric magnetic field (dipolar for example). It is quali- 
tatively clear that the slower the magnetic field decreases with 
distance, the more efficiently it accelerates the plasma due to 
a MC mechanism. The monopole-like magnetic field falls down 
while the dipolar magnetic field falls down faster than 
Correspondingly, the plasma is expected to be more effi- 
ciently accelerated in the monopole-like magnetic field in com- 
pare with the acceleration in the dipolar one. This is why the 
results obtained in the monopole-like model rule out efficient 
MC acceleration of the plasma in any other more realistic ax- 
isy mmetr ic model. Calculations performed by Contopoulos et 
al. (1999) confirm this conclusion. 

Thus, up to now all the attempts to propose an effective 
mechanism of relativistic plasma acceleration have not been 
successful. That is why the search for new possible mechanisms 
to explain relativistic plasma acceleration remains one of the 
most important problems of relativistic astrophysics. 

The disappointing conclusion regarding the low efficiency 
of the MC acceleration is valid only for axisymmetric magnetic 
fields, as previous studies considered this process only within 
the frameworks of axisymmetric models. However, the mag- 
netic field of real astrophysical objects is not axisymmetric. 
Radio pulsars give us the brightest example of this relation. 
Therefore, it is important to know how the non axisymmetry 
of the magnetic field affects the process of MC acceleration. In 
this paper we try to answer this question. 

To clarify the role of the nonaxisymmetry of the magnetic 
field for the MC acceleration of plasma, a model with the ini- 
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tially monopole-like magnetic field is used in this paper. This 
model is often used to investigate the processes of plasma coUi- 
mation and accelerati on in the rotat ing rn agnetosphe re (M ichel 



196C , Mest el fc S elley 



skin et al. Jl998 
I 



1970, Sakurai 



1985 



Bogov alov 



1999 



Bogovalov & Tsinganos 
The word "initially" implies that the mag- 



1992 



Be- 

Tsinganos & 



Bogovalov 

netic field of the non rotating star looks like the magnetic field 
of the "magnetic monopole". This model is very convenient 
since there are no closed field lines in the fiow. It remarkably 
simplifies the analysis of the plasma acceleration. 

The main goal of this paper is to answer the question: how 
does the nonaxisymmetry of the magnetic field of a star affect 
the process of plasma acceleration? An attempt to solve the 
problem of the structure of the magnetosphere of an oblique 
rotator with a dipole magnetic field in mass-less approximation 
has been done by Beskin et al. ( 1983 ) . The first step in the so- 
lution of this problem in the MHD approximation has been 
done in the work of Bogovalov (1999). In this work, the prob- 
lem of plasma flow in the magnetosphere of the oblique rotator 
with an initially split-monopole magnetic field was solved. The 
modulus of the magnetic field was axisymmetric and only the 
direction of the magnetic field varied with time and the az- 
imuthal angle. It was found that the acceleration remains non 
efficient in this flow as well. In the present paper, we are in- 
terested in the affect of the nonaxisymmetry of the modulus of 
the magnetic field on the process of plasma acceleration. 

Plasma fiow is described by a system of non linear equa- 
tions in partial derivatives. Usually these equations can be in- 
tegrated only numerically. There is one exception from this 
rule. Sometimes the problem can be solved self-consistently 
and analytically if we are interested in small corrections to 
the known solution. These corrections can arise if the known 
fiow is perturbed slightly by additional small forces. In partic- 
ular, slow rotation can be considered as a small perturbation 
of some known fiow from the non rotating star. This approach 
was firstly successfully applied for the numerical self-consistent 
solution of the fiow of t he so lar wind by Nurney and Suess 
(1975). Later Bogovalov (1992) has used this approach for the 
fully analytical, self-consistent solution of the problem of the 
cold plasma outflow, relativistic as well as nonrelativistic, from 
a slowly rotating star. In recent years, this method has also 
been used to so lve a range of problems by Beskin (see Be- 
skin & Okamoto 200C ) and references therein) . The idea of the 
present work is the following: if the nonaxisymmmetry of the 
magnetic field modifies the process of MC acceleration of the 
plasma, this modification can be seen at the level of the slow 
rotation of the object. Therefore, to answer the question: does 
the nonaxisymmetry of the magnetic field modify the process 
of acceleration of the plasma, we investigate the acceleration of 
plasma in the nonaxisymmetric magnetic field at slow rotation. 

This paper is organised as follows. In Sect. 2 we present 
ideal MHD equations defining the dynamics of cold relativistic 
plasma. A self-consistent analytical solution of the problem is 
given in Sec. 3. The solution in the wave zone is given in Sect. 
4. The possible astrophysical implications of the results are 
discussed in Sec. 5. 



2. Basic equations and relationships 

The system of time dependent equations defining the temporal 
evolution of the relativistic plasma outflow in an ideal MHD 
approximation is as follows (Akhiezer et al. 1975): 
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mn(^^ + (vV)7v) = g- E+ -!-jxH, 

at c 



(1) 



— — — —curl hi, 
cdt 




(2) 


, Tx 47r . 

curl H — — J + 
c 




(3) 


div H = 0, 




(4) 


E = 4nq, 




(5) 


dn ,. , , 
— — h divrvn) — 
at 


0. 


(6) 



E+ iv X H = 0. 



(7) 



Gravitation is neglected in these equations. It is convenient 
to consider the plasma flow in a spherical coordinate system. 
In this system, the equations of motion are 

mn(—^ + {v\/)'yvr — —) = 

1 r 1 dHr ld{rH^) 



rsin6' dip r dr ' 
djrHe) OH 1 dE^ dEe ^ 



(8) 



mn{ 



d^ve 



+ (vV)7ue + 



^{vrVe - cot 61) , 



1 ( d(rHo) dHr 



. d sin OH^p 



de 



dHe 



+ 



1 dEr dE^ I 



(9) 



mn{^^ + {vV)jv^ + — — - — ^ -] 



E^ 



47r L r sin 



1 ^dsinOHp dHg 



de 



dip 



1 dHr drHp ^ 
rsinO dip rdr ^ 

and the induction equations: 

dHr _ 1 a(sin6>g^) _ dEe , 
cdt rsinO rdO dp ' 

dHe _ 1 dEr _ d{rE^) 
cdt r sin 9 rdp rdr ' 



dH^ 
cdt 



, d{rEe) _ dE^, 
^ rdr rdO' 



(10) 

(11) 
(12) 
(13) 



The laws of conservation of the magnetic field and the mat- 
ter flux, together with the Coulomb law, take the following 
form: 



1 d{r^Hr) 1 d sin OHg 1 dH^ 



= 0. 



dr rsinO d6 rsin6' dip 

dn 1 dir^nvr) 1 9sin^nu9 1 dnv^ 



+ 



+ 



dt r^ dr r sin 9 d9 r sin 9 dp> 

1 d{r^Er) 1 d sin 9 Eg 1 dE, 



(14) 
0.(15) 



+ 



+ 



dr rsin9 d9 rsin9 dp 

The ideal MHD conditions are 



Er + -{vgH^ - VpHg) = 0, 



Eg + - {v^Hr - VrH^) 



0, 



E^ + -{VrHg - VgHr) = 0. 



= 47rg. (16) 

(17) 
(18) 
(19) 



In this system, q is the induced space electric charge den- 
sity, 9 is the polar angle and p is the azimuthal angle. 

The electric field in the axisymmetrically rotating steady 
state magnetosphere is connected with the poloidal magnetic 
field as follo ws: E = rsin^fi/cBp x and E^ = (Weber & 
Davis 1967). For the nonaxisymmetric magnetic field the same 
relationship is also valid, provided that the flow is in the steady 
state and all the variables vary periodically with time (Beskin 
et al. 1983). In the spherical coordinates, the relationships take 
the form 



rsin^Sl 

Eg = Hr 



Er = — '■ Hg, Etf, = 0. 



(20) 



It is easy to show that these relationships indeed satisfy all 
the equations of the electric field in a rotating magnetosphere. 
At stationary rotation, all the variables depend on ip and t in 
the spherical system of coordinates as the difference ^ = ip—Qt. 
Then, the induction equation ( |ll| ) is reduced to the equation 



9^ r sin 9 



+ 



= 



(21) 



which is fulfilled due to Eq. ( po| ) . The same is valid for the 
induction equation (|l2[) which is reduced to the equation 



d^ Er 



Q.Hg 



d^ r sin 9 



)=0. 



(22) 



The third induction equation ( |13[ ) after the substitution of 
Eqs. die] ) is reduced to the magnetic flux conservation condition 
(p^. Eq. (^) also satisfies the boundary conditions on the 
surface of the star where the tangential components of the 
electric field should be continuous. It follows from the ideal 
MHD conditions ( [l7| - p^ that below the surface of the star 
the tangential components of the electric field are E^p = 
and Eg — —Qrsin9/cHr. They exactly equal the tangential 
components of the electric field directly above the star surface, 
given by (^. 

3. Acceleration of the plasma at slow rotation 

3.1. Acceleration in the axisymmetric magnetic field 

The solution of the problem of steady state plasma outflow can 
be expanded in powers of SI. The poloidal magnetic field, the 



4 



plasma density, the poloidal velocity and the Lorentz factor of 
the plasma are invariant in relation to the rotation direction 
at the axisymmetric outflow. Therefore they are expanded in 
even powers of Q. Correspondingly, the toroidal magnetic field 
and the toroidal velocity are expanded in odd powers of Q, 
since these variables change sign at the change of the direction 
of rotation. Therefore, the expansion of the solution in the 
powers of Q takes the form 



Vp,o + Vp,2 + . 



n = no + fl 712 



7 = 70 + f2^72 + f2*74 + . 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



where lower index "0" denotes the variables for the non- 
rotating star. The first corrections to the unperturbed cold 
plasma outflow from the nonrotating star arising at slow ro- 
tations of the star with a monopole-li ke m agnetic fleld were 



calculated analytically by Bogovalov (1992). The first order 
corrections to the toroidal magnetic flel( 
locity are as follows 



d and the toroidal ve- 



rsinO 



where 



H 



p,0 



(29) 



(30) 



is the monopole-like magnetic fleld of the star and the lower 
index "*" denotes the values on the surface of the star. It fol- 
lows from (^9|) that the flrst correction to the energy of the 
plasma accelerated in the monopole-like magnetic field is pro- 
portional to fi^. Indeed, the equation for the Lore ntz factor of 
the plasma has the form (Landau & Lifshitz 1975) 



mnc'(§ + (vV)7) =jE. 
at 



(31) 



The second order correction to the Lorentz factor is defined 
as follows 



2 r,2d 
mnoc vqH —72 
r 



rQsme d . , 
An rdr 



(32) 



But according to ( p9| ) and ( po] ) the right hand part of this 
equation is equal to zero. Thus, 72 = 0. The fact that the first 
correction to the Lorentz factor at the axisymmetrical outflow 
begins with the term proportional to 0^ reflects the very low 
efficiency of the MC acceleration in the monopole-like magnetic 
field. 



3.2. Acceleration in the nonaxisymmetric magnetic field 

The Lorentz factor of the plasma accelerated in the nonax- 
isymmetric magnetic field is already not invariant in relation 
to the direction of the rotation. In this case it can (but not 
must) depend on the sign of fl. The expansion of the Lorentz 
factor in $7 takes the general form 



7 



: 70 4- r27i 4- n^72 + n^73 + ••■ 



(33) 



In the axisymmetric case, 71 = 73 = (generally speaking, 
the term 72 = applies only for the monopole-like magnetic 
field). In the nonaxisymmetric case these terms can (but not 
must) differ from zero. If for example 71 7^ 0, it would imply 
that the acceleration of the plasma in the nonaxisymmetric 
magnetic field is essentially more efficient in comparison to the 
acceleration in the axisymmetric magnetic field, at least for 
the case of slow rotation. To make sure that indeed 71 7^ 0, 
it is sufficient to do this for the magnetic field having small 
nonaxisymmetry. In this case, the problem can be solved ana- 
lytically. 

3.3. Outflow at fl = in the slightly nonaxisymmetric mag- 
netic field 

Let us initially consider the cold plasma outflow from the non- 
rotating star with a slightly nonaxisymmetric magnetic fleld. It 
follows from ( |3l| ) that the energy of the cold plasma is constant 
since the electric field E = everywhere. For the slightly non- 
axisymmetric magnetic fleld the solution also can be expanded 
on a small parameter, e, characterising the nonaxisymmetry of 
the initial split-monopole magnetic fleld. This expansion has 
the form 



Hp — Hp,o + ehp.i + 

Hip = shp^i + £^/itp.2 + ... 

Vp = Vp,o + eVp,i + e^Vp,2 

n = no + ehi + 

but 
7 = 70. 



(34) 
(35) 
(36) 
(37) 
(38) 

(39) 



since the energy of the cold plasma is constant. It is conve- 
nient to consider the perturbation of the magnetosphere due to 
a non uniform distribution of the magnetic field on the stellar 
surface. We assume for simplicity that the magnetic fleld has 
its own axis of symmetry m not coinciding with the axis of 
rotation, as is shown in Fig. ^. It is convenient firstly to solve 
the problem in the coordinate system where the magnetic fleld 
is axially symmetric and then to transform the solution to the 
laboratory coordinate system. The equation for the first cor- 
rections in the system where the magnetic fleld is axisymmetric 
follows from Eq. (^) and takes the form 



mnojovo- 



drVy 



rdr 



Ho ^ drh^^i 
4m rdr 



dhr,i 
rdx 



(40) 



where x is the polar angle in the system of the symmetry of 
the flow. It is convenient to introduce the flux function ■(/). The 
components of the poloidal magnetic field can be expressed 
through the flux function as follows 
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H 



Hr 



^ r sin X t^'" ' ' sin x "^^X 
The expansion of the function ^ in powers of e is 

= "Cosx + e/(r) + ....)• 



(41) 



(42) 



Taking into account that due to the ideal MHD condition 
V = vtl/H, it is easy to obtain the equation for / (Bogovalov 
1995) 



frr + fr 

r 



J.2 Q2j^ 



(43) 



where rj = cosx, Ra ~ \/ HfRi /47rmn*7oi'Q is the Alfven 
(or fast magnetosonic) radius of the unperturbed flow. The 
perturbation / can be presented as 



(44) 



where Qm are the eigenfunctions of the differential equa- 
tion (1 — Tf)d^ f /d^r)Qm ~ —m{m + l)Qm- They have the 
form Qmiv) = (1 ~ 'n^)9Pm{ri)/dVj where Pm are the Legendre 
polynomials of the order m. 

The equation for fm takes the form 



2x 



dfm 
dx 



m{m + 1) 



/.n = 0, 



(45) 



where x — r/Ra- The general solution can be presented as 
a sum fm = amZm{x) + bmYm{x) of two independent solutions. 
One of them, Zm, is regular at the point x = 1 and another, 
Ym, is singular at this point. The condition of regularity of the 
solution at this point gives bm = 0. The regular solution of Eq. 
( ^sl ) is Zm{x) = Pm{i/x). Therefore, the general solution of 
the problem regular at the Alfven surface can be presented as 
follows 



It is seen that the radial perturbation of the magnetic field 
is positive at the magnetic poles and negative at the equator 
(at e > 0). It is interesting that the perturbation of the radial 
component of the magnetic field decreases as at r — * oo 
. The initial magnetic field has the same dependence. This 
means that the fiow is not spherically symmetric at r — > oo. 
This is one of the features of supersonic fiows. The magnetic 
forces decrease with distance so fast that the supersonic plasma 
moves ballistically along straight lines defined by the initial 
conditions. The magnetic pressure gradient across the fiow lines 
is not able to change the trajectories of the particles of the 
plasma to restore the spherical symmetry at r ^ oo. 




X 



Fig. 1. The geometry of the coordinate systems connected with the 
axis of rotation (axis f2) and the axis of symmetry of the magnetic 
field (axis m) inclined to the rotation axis f2 under the angle a. 
is the unit vector directed along r. 



E 



{l/x)Qmiv)- 



(46) 



The boundary conditions on the stellar surface define the 
numerical coefficients am- 

It is convenient to take one term from expansion ( ^ ) which 
corresponds to the simplest non uniformity of the magnetic 
field on the stellar surface. It is especially interesting to con- 
sider the non uniformity which increases the magnetic field 
at the magnetic poles and decreases the magnetic field at the 
magnetic equator. The term with m = 2 corresponds to this 
distribution. The correction to the solution in this case is 



_ (3/3:^ - 1) . 2 

f = , , „ sm Y cos Y. 

(3/xi - 1) ^ 

The corrections to the magnetic field are 



hr,l 



6H,Ri 



Rlx^i/xi - 1) 

(3/x^ - l)H,Rl 
r2(3/x2 -1) 



sm X cos X, 



(2 cos X - sin x)- 



(47) 

(48) 
(49) 



The components of the arbitrary vector A in the spherical 
coordinate system with the polar axis directed along the sym- 
metry axis of the magnetic field m (dashed coordinate system) 
can be transformed into the components of the same vector in 
the laboratory coordinate system, shown in Fig. |l| as follows: 



Ae 



A^ (cos a sin — sin a cos 6 cos ip) — A^^, sin ip sin a 
sinx 



A^ (cos Q sin 9 — sin a cos 9 cos (f) + A^ sin (p sin a 
sinx 



where a is the inclination angle of the magnetic axis in relation 
to the axis ft. The equation connecting the coordinates in the 
systems is as follows 



cos X = sin a sin 6 cos ip + cos a cos t 



(50) 
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7r/2 (orthogonal rotator) is of special interest for us below. The 
components of the vector are connected as follows in this case 



Ag = -A^ 



and 



A^p — A^ 



cos t) COS ip 



\/l — sin^ 9 cos^ (f 



sm ip 



\/\ — sin^ Q cos^ LP 



(51) 



(52) 



3.^. Plasma acceleration in the slightly nonaxisymmetric mag- 
netic field at slow rotation 

At slow rotation and small nonaxisymmetry of the magnetic 
field the solution can be expanded in two parameters: Q, and 
e. The expansion takes a form 



Hp — Hp^o + ehp^i + .... 

= flHip^i + ehip.i + 
Vp = Vp,o + eVp,2 + 

n — no + eni + 

7 = 70 + eO.S'y + 



It follows from (^) that satisfies the equation 



2 ^^^7 
mnoc Qevo—— = Egjg, 
or 

where 

c 1 dhr.i d 

JO = £-ri — ^rhy,,i), 

47r rsmt? 0(p ror 



(53) 
(54) 
(55) 
(56) 
(57) 
(58) 

(59) 

(60) 



and Eg is defined by Eq. (|^). Let us calculate S'y for the 
particular case of orthogonal rotation (a — n/2) for the flow 
in the equatorial plane {6 = 7r/2). The transformation of the 
solution (pll p3) to the laboratory coordinate system gives 



c QeHt, sin ^ cos ^ 
4^ r^{3/xl - 1) ■ 

Integration of Eq. (|9|) with Eq. ( pl[ ) gives 

flR* vq 6esin5cos^ 



7 = 70 70 - 

where 

ip-Q.t. 



{3-xl) 



(61) 



(62) 



(63) 



It is seen that the violation of axial symmetry of the mag- 
netic field essentially modifies the acceleration process. The 
first correction to the Lorentz factor of the plasma is propor- 
tional to in the axisymmetric magnetic field, while even the 
small nonaxisymmetry of the magnetic field makes the depen- 
dence much stronger. The first correction to the Lorentz factor 
begins with the term proportional to Q. 



acceleration 




Fig. 2. Outflow in the magnetic field of a star with a monopole-Uke 
magnetic field having small non uniformity on the surface of the 
star. The magnetic field is slightly stronger at the poles and slightly 
weaker at the magnetic equator. The magnetic field is axially sym- 
metric in relation to the axis m. The particular case when the axis 
of the symmetry of the magnetic field is perpendicular to the axis 
of rotation (orthogonal rotator) is presented. At a slow rotation, the 
acceleration of plasma occurs along the field lines which curve with 
a radius Rc directed into the direction of rotation shown by thick 



The schematic structure of the field lines unperturbed by 
the rotation is shown in Fig. ^ The view from the top is pre- 
sented. The acceleration of the plasma takes place in the sectors 
where sin ip cos ip < (at t = 0, (y3 = at the direction along 
m). In other sectors the plasma is decelerated. It follows from 
the figure that the plasma is accelerated where the field lines 
unperturbed by the rotation are curved into the direction of 
the rotation and is decelerated where they are curved in the 
opposite direction. This behaviour of the plasma has a simple 
mechanical analogy. 

The motion of the magnetised plasma can be considered 
in somesense as a motion of a bead on a wire (Blandford & 
Payne 1982). The bead is accelerated centrifugally during the 
motion along the wire, provided that the wire is not in the 
shape of an Archimedian spiral with the components dl^p and 
dlr defined by the equation 



dlr 



It is easy to demonstrate that any departure of the wire shape 
from the Archimedian spiral results in variation of the energy 
of the bead. The bead will be accelerated centrifugally un- 
der the condition \-^\ < \ ''^J'" ^ | and decelerated in the op- 
posite case. The plasma behaves exactly in accordance with 
this mechanical model. It follows from (291) that the field lines 
of the axisymmetric magnetic field form the Archimedian spi- 
ral at small fl. At the nonaxisymmetric outflow, the toroidal 
magnetic field is a linear superposition of perturbations from 
rotation and the nonaxisymmetry. The initial nonaxisymme- 
try straightens the magnetic field lines in the sectors where 
simp cos ip < 0. The shape of the field lines here satisfies the 



inequality 



< 



The plasma is accelerated in this 



case. In other sectors, the field lines are curved in an azimuthal 
direction more strongly in comparison to the Archimedian spi- 
ral. Therefore, the plasma is decelerated. 
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4. Generation of magneto sound waves 

The modification of tlie MC acceieration is not tfie only new 
feature of the outflows in the magnetic field with violated 
axial symmetry. The source with the nonaxisymmetric mag- 
netic field excites the magnetohydrodynamical (but not mag- 
netodipole) waves in the wind. Perturbation theory allows us 
to investigate this process while the amplitude of these waves 
is small. 

In the ajdsymmyteric case, the azimuthal magnetic field 
generated by the rotation is defined by Eq.(|29|). Perturbation 
theory can be applied while the pressure of the toroidal mag- 
netic field in the co-moving coordinate system is less than the 
press ure of the unperturbed poloidal magnetic field (Bogovalov 
1997). In the laboratory coordinate system this condition takes 



the form 



< Hp, 



(64) 



This inequality limits the application of the perturbation 
theory by the distances r less than ^qVo/^- When the last 
value becomes less than Ra, the perturbation becomes strong 
everywhere in the flow. Thus, the perturbation theory can be 
applied only if 



70 I'D 



< 1. 



(65) 



The objects satisfying condition ( |65| ) are referred to as slow 
rotators. In the opposite case they are fast rotators. Analysis 
shows that all radio pulsars, including the Crab pulsar, are 
slow rotators at the conventional parameters of the winds (Bo- 
govalov h999|) . 

Eq. (|6E|) defines the condition under which the perturbation 
of the flow in the subsonic region remains small compared with 
the initial magnetic field. Our use of the perturbation theory is 
valid given an additional restriction. During the calculation of 
the first corrections, we neglected the contribution of the en- 
ergy density of the toroidal and electric fields generated by the 
rotation into the inertial mass of the plasma. In order words, 
it was assumed that the condition 



> 



An 



(66) 



is fulfilled, where the upper index " " denotes values in 
the co-moving coordinate frame. It results in the additional 
condition 



Ra «C Rl, 



(67) 



where Rl = i^/c is the radius of the light cylinder. Thus, 
our consideration is valid under the conditions (^) and ( p7[ ) 
in the region limited by the radius 7o«o/f^- 

The objects with weak violation of the axial symmetry ex- 
cite small amplitude MHD waves in the outfiowing wind, which 
can be considered in the framework of the perturbation the- 
ory. These waves are formed in the wave zone with dimen- 
sions of the wavelength A = Vm/fi (Landau & Lifshitz 1986), 
where Vm is the velocity of the wave propagating in the mov- 
ing plasma. The fast magnetosonic velocity of the cold wind 
decreases with an increase in r. In the region r > Ra, the fast 
magnetosonic velocity of the plasma Vf is small compared with 



«o for slow rotators and the condition vo ^ Vf is fulfilled ev- 
erywhere along the flow, under condition ( |67| ) . The wavelength 
can be estimated as A ~ Vo/Q. in this region. Thus, four scales 
appear in the problem. One is the initial radius of the Alfven 
(or fast magnetosonic) surface, Ra- Others are the radius of 
the wave zone A, the radius of the light cylinder Rl and the 
radius Vmax = 7oA, where the perturbation theory itself can 
not be applied. The inequality -Ra/7oA < 1 is fulfilled for the 
slow rotators. For the nonrelativistic plasma, Vmax = A. The 
wave zone is located in the region where perturbation theory 
can not be applied. Therefore, the generation of MHD waves 
in the nonrelativistic winds cannot be considered in the frame- 
work of the perturbation theory developed here. However, in 
the relativistic plasma with 70 ^ 1 there is a gap between 
A and 70A where the perturbation theory remains valid. The 
MHD wave generation can be considered in the framework of 
this theory. 

It is easy to extend our s olutio n to the zone A < r < 70A, 
following Landau & Lifshitz (1986, sect. 74). The electric and 
magnetic fields in the outgoing wave depend on r and t as 
Ai{t — r/vo)/r, where Ai are unknown functions. In the zone 
R, <^ r < X the solution has a form Bi{t)/r. Both solutions 
must coincide at r = A. This implies that Ai — Bi(t — r / vq) . 
So, we need only to define functions Bi of the solution obtained 
in the zone r < A to obtain the solution at r > A. 

It was pointed out that the perturbation of the radial com- 
ponent of the magnetic field hr, proportional to e, falls down 
as at large r. This perturbation results in the perturbation 
of the toroidal magnetic field (see Eq. (29)) 



Q,r sin 6 



(68) 



and in the perturbation of the 6 component of the electrii 
field (see Eq. (20)) 



Q,r sin 9 , 



eg 



(69) 



The electric and magnetic fields from ()^^ and (^9|) decrease 
with distance as at r ^ 7?«, since hr has the following 
dependence ( see Eq. (49)) 



hr 



H,Rl 



>2(3/a;2 



-(3cos X - !)• 



(70) 



All other terms in Eq. (|4£|) are neglected since they decrease 
with r faster than . According to Eq. ( po[ ) and Eq. (|6^), 
cosx depends on t and tp in the combination £^ = ip ~ Q,t. To 
obtain the solution at r > A, we have to replace it by the 
combination tp — ^{t — t/vo). As a result, cosx takes the form 
(see Eq. (g)) 



cos X = sin a sin 9 cos{lp — Q.{t — r/vo)) + cos ( 



(71) 



The solution ( |68| , p9y7C| , |7lD allows us to define a correction 
to the rotational losses of the central object due to the emission 
of the MHD waves in the wind and to estimate the effect of 
the wind acceleration by these waves. The correction to the 
r component of the Poynting flux S = -^E x H due to the 
perturbation of the electric and magnetic fields is as follows 



Sr = —evQ 



HlRl 
27rr2(3/a;2 



,^R,sin9 2,„ 2 
-) (3 cos X 



"0 



(72) 
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The ratio of this correction over the particle flux density in — _5£( 



mc units is as follows 



i?aJ^ sin S ,2 sin 6 cos X sin a sin 5 rO, 

c ' (S/xl - 1) 



(77) 



5a = -2e( 



QRaam9.2 (3cos^x~l) 



70 



{3/xi 



(73) 



This value shows the maximal possible variation of the 
Lorentz factor of the plasma at the full absorption of the MHD 
wave. 

Integration of the correction to the Poynting flux over a 
surface surrounding the source gives the first correction to the 
rotational losses of the central object, 



SE - ^ ,n^^El^^ 



(3 cos a — 1). 



(74) 



The zero order expression for the rotational losses is given 
in (Bogovalov 1999). As was expected, the departure from uni- 



form distribution of the magnetic field lines on the surface of 
the star changes the rotational losses. Correction ( |74| depends 
on the sign of e and does not go to zero even at axisymmet- 
ric rotation (a = 0). This fact has a simple explanation. The 
total Poynting fiux from an axisymmetrically rotating object 
depends on the distribution of the poloidal magnetic fi eld o n 
latitude at a large distance from the source (Bogovalov 1999). 
At the fixed total magnetic flux, the rotational losses increase 
when the magnetic flux is more concentrated at the equator. 
The concentration of the magnetic fiux at the poles conse- 
quently results in the decrease of rotational loss. Eq. (Q) de- 
scribes this dependence. The magnetic fiux is concentrated at 
the poles for e > on the surface of the star, but according 
to (49), at large distances the magnetic fiux is more concen- 
trated at the equator. Therefore, the rotational losses increase 
for e > and decrease in the opposite case. 

Qualitatively, this also explains the fact that according to 
Eq. ( |74[ ) the rotational losses decrease at the orthogonal rota- 
tion (at e > 0). In this case, the magnetic fiux at large distances 
is more concentrated along the plane of the magnetic equator, 
oriented perpendicular to the rotational equator. The average 
magnetic field in the equatorial plane decreases and results in 
the decrease in rotational losses. 

The MHD wave propagating in the wind changes the energy 
of the particles of the plasma. The variation in the energy is 
described by Eq. (31) with the electric current j = -^[curlil — 

- ) . Substitution of solutions (psl Fil) in this expression gives 



Je = 



4-K 



2 2 

«o7o 



(75) 



The correction to the Lorentz factor of the plasma depends 
on its arguments, as Sj{r,^). After substitution of this expres- 
sion and ( [zil) in Eq. (31), we obtain 



85"/ 
dr 



-6eR 



2 SI sin d cos X sin a sin ^ 



-yoVQC^ 



{■i/x. 



(76) 



Formally, integration of Eq. ([76|) over r gives an infinite 
result. It is clear why this should occur. We do not take into 
account the decrease in the amplitude of the MHD wave due to 
the redistribution of energy between waves and particles of the 
plasma. However, we can estimate the variation in the energy 
of the plasma depending on r. Integration of Eq. ( fz^ ) gives 



The lower limit of integration is neglected here. This cor- 
rection can be negative or positive, depending on the sign of 
cos X sin ^. Thus, we consider the process of redistribution of 
ener gy, but not pure acceleration. Comparison of Eq. (^^ with 
Eq. ( |73[ ) shows that even at the distance 70A only a small part 
of the energy of the wave is transformed into the energy of the 
plasma and vice versa at 70 ^ 1. So, this process is extremely 
noneffective in the relativistic regime. 

5. Implications for the physics of relativistic winds 
from radio pulsars 

It is generally believed that plasma is produced and ini- 
tially accelerated in the pulsar magnetosph ere in the so- called 
electrostatic gaps (Ruder man fc Sutherland 1975, Arons 1981, 
Cheng et al. |l984 Romani p9|). The ideal MHD does not ap- 



ply to these gaps. However, these gaps occupy only a small part 
of the magnetosphere. They produce the initially relativistic 
plasma, which is dense enough to screen the electric field and 
to provide the ideal MHD conditions outside the gaps. There- 
fore, the fiow of this plasma can be described in an ideal MHD 
approximation as the wind. Below, we assume for simplicity 
that the wind with a prescribed initial density and Lorentz 
factor is formed somewhere near the surface of the pulsar. 




Fig. 3. The view from the top of the schematic magnetic field struc- 
ture of the outflow from an orthogonal rotator with an initially dipole 
magnetic field. The light cylinder is shown by the solid circle. Di- 
rection of rotation is shown by the thick arrow. The region of the 
closed field lines does not reach the light cylinder and is limited by 
the dashed circle. Some of the open field lines (shadowed) are curved 
in the direction of rotation up to the light cylinder. The magnetocen- 
trifugal acceleration will take place on these field lines. Beyond the 
light cylinder, all the open field lines are curved in the Archimedian 
spiral typical of axisymmetric outflows. 

The schematic structure of an orthogonal rotator with an 
initially dipole magnetic field is presented in Fig. ^ The dipole 
magnetic field is not axisymmetric at oblique rotations. There 
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are regions of closed field lines in the magnetosphere and re- 
gions of open field lines along which the plasma outflows. Some 
of these field lines are initially curved in the direction of ro- 
tation. Therefore, the new type of MC plasma acceleration 
discussed in this paper must work under radio pulsar condi- 
tions. The question remains, however, whether this mechanism 
of plasma acceleration essentially changes the energy of the 
plasma in pulsar conditions. 

Real pulsar parameters corresponds to the conditions Ra > 
Rl and Ra < JoRl- So, our solution cannot be directly ap- 
plied to radio pulsars since it is valid at Ra < Rl- However, 
this solution can give important qualitative information about 
the process of wind acceleration in any magnetosphere having 
violated the condition of axial symmetry. 

To understand at least qualitatively the influence of the 
violation of the axial symmetry on the process of the MC ac- 
celeration in a pulsar magnetosphere we should estimate the 
dimensionless values of the expansion parameters in Eq. ( p2| ) 
for pulsar conditions. If these parameters are small, it would 
imply that this effect is not important for radio pulsars. If these 
parameters are of the order of 1, it would imply that the MC 
acceleration in the pulsar conditions can essentially change the 
energy of the plasma. 

The expansion in Eq. (^2|) occurs in two parameters e and 
^^S^. £ can be estimated as e ~ SHt/ < Ht > where SHt is the 
maximal variation of the magnetic field on the surface of the 
star and < iif* > is the average magnetic field on the surface 
of the star. It is evident that e ~ 1 for dipolar magnetic fields 
which are believed to occur on the pulsar surface. 

In the second dimensionless parameter the physical 

sense of R, is the most unclear. At first glance, it seems that 
this parameter defines the scale over which the magnetic field 
decreases with distance from the star. However, more close in- 
spection of Eq. (^2|) shows that this cannot be correct. The first 
correction to the Lorentz factor (|6^) does not depend on the 
magnetic field _ff* at all, provided that 3Ra ^ Rl- Therefore, 
_R* cannot be interpreted as describing the scale of decrease of 
the magnetic field. 

To understand the physical sense of _R* in Eq. (^), it is 
necessary to note that the MC acceleration under considera- 
tion occurs only in regions where the magnetic field lines were 
initially curved into the direction of rotation. Simple analy- 
sis of the correction to the mag netic field @ shows that the 
violation of axial symmetry on the surface of the star results 
in the bending of field lines in the region with a character- 
istic size R,. The shape of the field lines does not depends 
on H, at 3Ra ^ R^- This inequality means that the energy 
density of the magnetic field exceeds the energy density of the 
plasma near the star surface. At larger distances, the field lines 
straighten and do not accelerate the plasma. Thus, R, plays 
the role of the characteristic scale Rcurve on which the mag- 
netic field lines have nonzero curvature. Only accidentally does 
Rcurve = R* in our specific case. Thus, the second parameter 
defining the expansion of the terminating Lorentz factor can 
be presented as ilRcurve/c. 

In the pulsar magnetosphere, the scale of the poloidal mag- 
netic field variation is equal to the radius of the star, _R* . The 
scale Rcurve depends on the topology of the poloidal magnetic 
field. In the conventional models, the region of the closed field 
lines almost reaches the light cylinder. Direct calculations in 
the a xisvrn metric model confirm this assumption (Contoupolos 
et al. 199£). For this magnetic field topology Rcurve ^ Rl and 



therefore, the parameter ilRcurve/c also should be close to 1. 
For this reason, we conclude that the influence of the violation 
of the magnetic field axial symmetry on the process of the MC 
acceleration is certainly not small in the pulsar conditions and 
it could provide efficient acceleration of the relativistic winds. 

It is difficult to estimate reliably the terminating Lorentz 
factor of the plasma accelerated due to this mechanism in the 
magnetosphere of pulsars. Any field line initially curved into 
the direction of rotation is transformed into an Archimedian 
spiral well beyond the light cylinder where acceleration does 
not occur. The maximal Lorentz factor is defined by the dis- 
tance where this transition occurs. This is defined by the field 
structure at the vicinity of the light cylinder near the last closed 
field line of the magnetosphere, as is shown in Fig. ^ (shadowed 
region). Only accurate mathematical solutions will give an ac- 
curate answer to questions regarding the terminating Lorentz 
factor of the plasma. However, it is easy to understand that if 
this mechanism really provides transformation of the Foynting 
flux into kinetic energy of the plasma, then the acceleration 
must take place very close to the light cylinder. 

In conclusion, we note that the acceleration of the relativis- 
tic winds at the light cylinder of radio pulsars looks attractive 
from an astrophysical point of view as well. The mechanism of 
MC acceleration of the plasma which we have described can 
allow us to resolve a long-standing problem of additional ac- 
celeration of the wind from the Crab pulsar. It follows from 
the comparison of the observations of the Crab Nebula with 
the MHD theory of the interaction of the pulsar wind with the 
interstellar medium that the ratio of the Foynting flux to the 
density of the kinetic energy fiux in the wind bef ore th e termi- 
nating sho ck is o f the order ~ 10~^ (Rees & Gunn 1974, Kennel 
& Coroniti|l98^). However, this ratio cannot be achieved at the 



dissipativeless axisymmetric outflow(see in particular Chiueh 
et al.'s (1998) discussion of this question). It was argued re- 
cently by Begelman (1998) that the ratio of the Foynting fiux 
to the kinetic energy fiux of the order of 1 can also be consis- 
tent with the observations of the Crab Nebula if we take into 
account a possible instability of the fiow after the terminating 
shock. In any case, an additional acceleration of the wind of 
the relativistic plasma is needed to explain the observations of 
the Crab Nebula, since the electrostatic gap models are not 
able to repr oduce a wind with the necessary characteristics 
(Arons 1996). If we assume that the wind is accelerated due to 
the MC mechanism and that it carries off approximately half 
of the total rotational energy of the pulsar in the form of the 
kinetic energy, this resolves the problem of wind acceleration. 
The ratio of the Foynting fiux to the kinetic energy flux of 
the plasma appears to be of the order of 1 immediately after 
the light cylinder, which is consistent with the observations 
of the Crab Nebula according to Begelman's scenario (1998). 
This also resolves the difficulties with the mechanism of the 
wind acceleration at large distances from the pulsar proposed 
by Co roniti (1990) and recently revisited by Lyubarskii & Kirk 
(200C). There is simply no need to invoke additional mecha- 
nisms in this acceleration. 

Analysis performed by Aharonian & Bogovalov (1999) (see 
also Bogovalov & Aharonian (2000)) shows that the wind from 
the Crab pulsar with energetics comparable to the total rota- 
tional losses of this pulsar and accelerated near the light cylin- 
der should have a total particle fiux of ~ 10'*''part./s and an 
average Lorentz factor of ^ 10^ to be consistent with the lim- 
itations imposed by observations of the Crab Nebula in VHE 
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gamma-rays. It is interesting that the wind with these parame- 
ters indeed provides the total injection rate into the Crab Neb- 
ula necessary to explain the electromagnetic e mission in the 
wide range from radio to X-rays (Rees & Gunn L974). Thus, 
the idea that the mechanism of magnetocentrifugal accelera- 
tion of the plasma in the axially non uniform magnetic field 
can be responsible for the acceleration of the wind in the mag- 
netosphere of radio pulsars does not contradict observations of 
the Crab pulsar and Nebula and deserves further development 
in a more realistic model. 
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